Field-induced quantum soliton lattice in a frustrated two-leg spin-Y2 ladder 
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The field-induced quantum phase transitions (QPT) of the spin ladder material B^Cui-^Zn^- 
POe have been investigated via 31 P nuclear magnetic resonance (NMR) on single-crystal samples 
with x = and x = 0.01. Measurements at temperatures between 0.25 K and 20 K in magnetic fields 
up to 31 T served to establish the nature of the various phases. In BiCu2P06, an incommensurate 
(IC) magnetic order develops above a critical field po-Hci — 21 T; the field and temperature depen- 
dences of the NMR lines and the resulting model for the spin structure are discussed. Supported 
by results of Density-Matrix Renormalization Group (DMRG) calculations it is argued that the 
observed field-induced IC order involves the formation of a magnetic-soliton lattice. An additional 
QPT is predicted to occur at H > H c i- For x = 0.01, this IC order is found to be stable against 
site disorder, although with a renormalized critical field. 

PACS numbers: 75.10.Pq, 76.60.-k, 05.10.Cc, 75.40.Cx 



Quantum domain walls (QDW) are an important con- 
cept for understanding the magnetism of low-dimensional 
systems [![. In one dimension (ID), QDWs are due to the 
fractionalization of a spin-flip excitation into two spin- 
Y2 objects and represent kinks separating two energet- 
ically degenerate domains. Since their description in- 
volves non-linear differential equations [H , they are usu- 
ally termed solitons, in analogy to those known in con- 
densed matter in relation with charge- and spin-density 
waves d, H|, lattice dislocations [3 or superfluid 3 He 
Q. In ID spin-!/2 antifcrromagnetic (AFM) Heisenberg 
chains with strong geometrical frustration or spin-phonon 
interactions, QDWs are non-zero energy excitations that 
divide the chain into two regions with different dimcr- 
ization patterns @~E1- Interchain interactions can select 
a dimerization pattern Q, whose dimer bonds can be 
broken by applying a magnetic field exceeding a critical 
value H c . A regular lattice of static QDWs then ap- 
pears. The interest in non-linear phenomena and the 
related non-trivial excitations motivates the search for 
solitons induced solely by str ong magnetic frustration, a 
natural cause of dimer order [1CJ, 1 . Field- induced (FI) 
QDWs were characterized in detail via NMR and neutron 
diffraction in the prototype spin-Peierls system CuGc03 
[12, . The scarcity of model materials in which dimer- 
ization is triggered exclusively by magnetic frustration 
has so far limited the number of related experimental 
and, to some extent, theoretical studies. 

Based on results of experiments and model calculations 



presented below, we argue that the spin- 1 /^ zigzag ladder 
compound BiC^POg, with geometric frustration caused 
by the next-nearest neighbor (NNN) interactions along 
the ladder legs, exhibits FI incommensurate (IC) mag- 
netic order involving the formation of a QDW lattice. 
Upon replacing Cu by Zn the ordered phase persists, but 
is distinctly affected by moderate quenched disorder. 

In BiCu2P06, the ladders formed by two coupled NNN 
frustrated chains run along the 6-axis of the crystal struc- 
ture ll5| . The inequivalence of the Cu 2+ sites allows 



the NNN exchange constants to alternate in magnitude 
[l| ( J2 and J2 in Fig.[IJ. Residual interladder interactions 
are expected along the c-direction, while PO4 tetrahedra 
act as non-magnetic spacers between different 6c-planes. 
The strong spin frustration 17[ and the unexpected ob- 
servation of a complicated set of first- and second-order 
FI transitions above the closure of the spin gap A [l|, [l8[ 
are reasons for the growing interest in this magnetic in- 
sulator. 

In the first part of this Letter we present 31 P NMR 
data obtained from 2x2x2 mm 3 single crystals [l5[ 
of B^Cui-^Zn^POg in magnetic fields up to 31 T, a 
range currently not accessible to, e.g., neutron diffrac- 
tion experiments. The field was applied along the &-axis 
with ±3° uncertainty and the covered temperature range 
was between 0.25 K and 20 K. The field dependence of 
the 31 P NMR lines in B^Cui-zZn^POe for x = and 
x = 0.01 is shown in Fig. 41] . We first consider the 
case of BiCu2P06- For H < H c \ the system is disor- 
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dered due to ID quantum fluctuations and the regular 
local magnetization results in the same sharp resonance 
for all the 31 P sites. For if > H c \ instead, a triple-peak 
NMR lineshape appears. It indicates a continuous distri- 
bution of local fields, suggesting magnetic order with an 
incommensurate spin pattern [4j. In general, the Cu 2+ 
electronic spin S^j contributes to the local-field experi- 
enced by the NMR nuclei via a priori unknown hypcrfinc 
interactions. In our case, a reliable evaluation of all these 
couplings was obtained by fitting the rotation patterns 
of the maxima of the main and satellite 31 P NMR lines 

for x = o.oi mim. Quantitative comparisons of spin- 
structure models with experimental data arc thus possi- 
ble. The four 31 P sites present in a unit cell are labelled 
by the index p. For H || b, the local field due to the i-th 
Cu 2+ ion at the j-th 31 P along the b direction can, ne- 
glecting a minor dipolar contribution, be approximated 

by h P U) * J2i,x AHS? + x,i - + £ ; A\{S\ +U + Si), 
where A = a, c, the ladder leg is denoted by / and the i 
values are integer multiples of 1 /2 (see Fig. [T] for the nota- 
tion) . The hyperfine couplings change sign according 
to symmetry when varying the index p. The longitudinal 
coupling A\ ~ 0.2 T/^ B is dominant. Since A" /c < A\ 
, the NMR spectra reflect mostly the longitudinal spin 
component S\ The normalized NMR lineshape 1(f) is 
dictated by the distribution of Jph p (j), with 7p if the 31 P 
Larmor frequency. The onset of FI order at H c i is mani- 
fest through the strong increase of both the first moment 
Mi and W, the square root of the second moment of 1(f) 
with field (see Fig. Hd) (H. For if < H cl both Mi and 
W show a linear increase with if, suggesting that S\i and 
5, ( A ; 7^ 0. In BiCu2PC>6, the 6-axis is a principal axis for 
the g-tensor and a possible Dzyaloshinskii-Moriya (DM) 
interaction associated with the Jj_ coupling must be par- 
allel to b If only these two terms are considered, the 
U(l) symmetry [2l| is preserved for H || b, leading to 
Mi = W = for if < H cl [H]. The fact that this docs 
not occur means that a DM interaction D associated with 
the J\ , J2 or J2 exchange paths is the leading anisotropic 
term. 

The onset of magnetic order is also indicated by dis- 
tinct maxima in the nuclear relaxation rate T^iT^H) 
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FIG. 1: Structural unit of a spin ladder in BiCu2P06, includ- 
ing the positions of 31 P nuclei. Only Phosphorus sites with 
p = 1, 2 are displayed. See text for details. 
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FIG. 2: (Color online) a) NMR lines /(/) in the field-induced 
phase of Bi(Cui_ I Zn a: )2P06. Lines are normalized to their 
peak maxima and their position is relative to the Larmor 
frequency. Data for both x = (foreground) and x — 0.01 
(background) at 0.25 and 1.4 K, respectively, are displayed. 
Magnetic field values, reported to the left (right) of the line 
for the x — (0.01) case, are relative to the critical fields 
[20.96(7) T and 24.21(9) T for x = and 0.01]. Arrows 
indicate the shift of the three blurred peaks for x = 0.01. b) 
Field dependences of the first (Mi) and the square root of 
the second moment (W) as computed from the lines in (a) 
for the x — case. Blue curves are guides to the eye, while 
ci = 4.1 x 10~ 3 MHz/T and c 2 = 0.73 x 10~ 3 MHz/T. 



[19( , whose (T, H )-dependences are displayed as red 
squares in the phase diagram of Fig. [3^,. The correspond- 
ing boundary was fitted using Tn oc [H c (T) — iJ c i(0)]*, 
resulting in <j> = 0.42(5); at T = K the onset of the FI 
order is at a critical field of fj, H cl (0) = 20.96(7) T. The 
value of the exponent 4> is in good agreement with a re- 
cent mapping of the phase diagram using torque magne- 
tometry at T < 1 K [23$ . Departures of 4> from (f> = 2/d, 
expected for rf-diniensional Bose-Einstein condensation 
of triplet excitations in FI order of quantum AFMs, may 
be caused by both anisotropy and frustration 2l|, |24 1 . 

The random suppression of spin degrees of freedom 
(known as site-disorder) via the introduction of spin- 
zero Zn for Cu at a few-percent level has a significant 
impact on all the regions of the BiCu2POg phase dia- 
gram, namely, the quantum-disordered region [k^T < 
A (if)], the IC phase [if > H cl (T)], and the quantum- 
critical point [H c i(x,T = 0)]. As shown in Fig. for 
x = 0.01 the critical field increases by about 3.2 T to 
fM>H cl (x = 0.01,0) = 24.21(9) T. Contrary to the case 
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of bond-disorder [25|, |27|, l28( , site-disorder allows for di- 
rectly probing finite-size effects. For B^Cui-zZn^POe 
a larger spin gap is the result of the finite-size ladders 
between the Zn sites, hence implying a higher critical 
field -f/ci(0.01, 0). FI magnetic order is absent in other 
disordered quasi- ID spin system with incommensurate 
correlations [25|, (2(|. Here, its existence can be demon- 
strated via NMR by monitoring the width w(T, H) of the 
NMR lines at 1/4 of the line maxima, which we choose as 
a measure of the order parameter. In Fig. [3]; we compare 
data obtained for H — H c i ~ and 2 T, respectively. The 
field-induced contribution is well visible above the Curie- 
Weiss type broadening. A triple-peak line shape, albeit 
blurred due to disorder, is also observed for x = 0.01 
(Fig. [2k). It indicates that the ladder regions between 
two Zn-sites tend to develop the same kind of IC order 
as for x = 0. Supporting this conclusion is also the ex- 
ponent <j> which, in the measured temperature range and 
within experimental error, is found to be the same for 
both x = and x = 0.01. 

Regarding the impurity-induced phase marked by dia- 
monds in Fig. [3^, we refer to T 1 _1 (T, H) shown in Fig.|3]3. 
For x = 0, an activated behavior reflects the gap in the 
excitation spectrum. Below T = 10 K there is a distinct 
upturn of Tj - , which is suppressed by the applied field. 
This can be related to residual impurity effects, as those 
occurring in gapped Haldane chains [29(. For x = 0.01, 
with controlled defects in the material, the same upturn 
is more pronounced. Here we observe maxima of these 
spin fluctuations occurring at a field-independent tem- 
perature, that disappear above 20 T. Most likely, the 
high-field IC-ordered phase is not directly adjacent to the 
low-field impurity-induced region, contrary to other Zn- 
doped cuprates possessing a spin gap [3(]| . This opens the 
possibility to experimentally investigate the new plateau 
phases proposed for site disorder [311 ] . 

In what follows we use numerical tools to study the 
IC spin structure realized in BiC^POg. We investigate, 
for H > H c i, the magnetization process of the ladder 
model with NNN exchange interactions proposed for Bi- 
CU2PO6 [H and show that it should not be seen as a 
simple increase of the number of individual bosonic-likc 
spin-1 triplon excitations, but rather as the formation of 
quantum domain walls (QDWs) with effective spin 1 /2. 
With our model, we can account for an additional QPT 
at H > H c i recently observed in experiments fl8| . 

From a classical point of view, the natural state of 
frustrated ID materials in a magnetic field is a canted 
spin helix 0,113. In a quantum-mechanical treat- 
ment, a renormalization of both the pitch angle and 
the magnitude of the ordered moments is expected 
[33| . Since in BiC^POg the geometrical frustration 
acts only along the ladder, we considered a tilted cone 
arrangement of the magnetic moments with Sn = 
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FIG. 3: (Color online) a) Phase boundaries of BiCu2P06 and 
Bi(Cui-. I Zn a: )2P06 for H || b. Points are determined from 
T^ 1 (T, H) relaxation data (see text for details), while <f> = 
0.42(5) and (pi = 0.41(3). b) T-dependence of Tf 1 for H < 
H c i, for x — (dots) and x = 0.01 (diamonds). The solid 
green line is a fit to cexp(-A/T), with A = 51.6(7) K. c) In- 
dependence of the NMR line widths w(T, H) at 1/4 of the line 
maxima for x = and x — 0.01. The expression c/(T + Tcw) 
represents a Curie- Weiss law, while the power-law exponent 
of the order parameter [fitted by fixing Tn from panel (a)] was 
found to be /? = 0.31(6) for the x = case. 



b by a polar and an azimuthal angle (77, e), while f3; = qi, 
with q || b as the propagation vector. For the or- 
dered magnetic moments m±, perpendicular to the cone 
axis, we used m\,± = — TO2,_i_. The resulting NMR 
lines, simulated for several 77, e, and q values, using the 
known hyperfine parameters, exhibit profiles as those in 
Fig. Hk, in clear disagreement with experimental data. 
The experimental resolution was taken as the width of 
the BiCu2P06 line at 19 T; the calculated line shift is 
/ = 2^ P {A\ + A^)m\i cos?y. 

In order to retain quantum effects in the model, we 
performed a DMRG analysis of the FI phase transi- 
tion in a system modeled by the Hamiltonian of two 
coupled frustrated J1-J2 chains, previously suggested 
for BiCu 2 P0 6 in Ref. Q (with J x = J 2 , J 2 = J 2 /2, 
J± = 3Ji/4, Ji/k B ~ 140 K). We consider the low 
triplet-density limit, i.e., H > H c i, in systems with 
sizes up to 2 x 128. The complete study will be pub- 
lished elsewhere 



34] 



tates the spin structure away from the applicd-ficld axis 



The main results reported here 
turn out to be qualitatively different from those ex- 
pected in the classical case. In a purely ID system of 
length L, one can draw a convenient analogy between 
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single triplet excitations and non-interacting hard-core 
bosons [35j. Therefore, with n = S^ Qt triplets in the 
finite system with open boundaries, the local magnetiza- 
tion is expected to be Sf iee l (i P ) = A liP J2k=i IV^p)! 2 
[3r3 |. where ipk{ip) is the particle-in-a-box wavefunction, 
as shown for the S* ot = 8 case in the upper panel of 
Fig. |4p. P is the site parity, i.e., the ip values are even 
or odd multiples of 1/2 and A ljevon/odd = A 2 , dd/even- Wc 
call such a high-field phase commensurate (C) because 
Sfree i(ip) and therefore mjn is constant in the system if 
L — > 00 and because, from numerical considerations [371 ] . 
we expect in that regime a spin helix with m>i t ± 7^ 
and a propagation vector close to the classical value 
q c = 7r- 1 arccos[-Ji/2(J 2 -|- J 2 )] = 0.608 « 2/3 0. The 
structure for q = 2/3 is sketched in Fig. 0b. Approach- 
ing the critical field from above, n decreases together 
with the overlap between particles, and the DMRG re- 
sults predict a C-IC transition at an L-dependent value 
of n, corresponding to a second critical field H = H C 2- A 
continuous QPT is expected in case of weak interladder 
interactions [34L l37|. 

In the IC phase, a triplet fractionalizes into two spin- 
1/2 objects. As shown in Rcf. 3J,|37J, these QDWs break 
the short-range dimer order present along the ladder legs 
in the large J± limit Q , instead of the better known long- 
range dimer order studied at the Majumdar-Ghosh point 
of J 1- Jochains for instance By analogy with the J± = 
case [Tj , we fit the S* ot = 1 data of Fig. 2b with the ex- 
pression S^jiip) = [1 + Bi :P sm(2niq + 9^p)}S^. cc l {ip), 
where q is incommensurate, 9^ CVC n = #z,odd + tt and 



ti/odd 



2, odd/even • 



In this model, the region 



of the phase diagram with FI order and H || b inves- 
tigated via NMR in the present work, corresponds to 
H c i < H < H C 2- For the comparison with experimental 
data we need to describe this intermediate phase in the 
L — > 00 limit. The latter is formed by commensurate re- 
gions between QDWs arranged in a regular pattern. For 
H > H C 2, contrary to the IC region 3, 3?}, the struc- 
ture is a canted helix such as the one shown in Fig. 0b. 
Hence we choose the equation for S,;.z to represent the 
untilted cone (77, e = 0) already used for the simulation 
in the left part of Fig. |4jt and use Bi = am[i/(k£), k] 
for the phase modulating the transverse component (see 
also [371]), In this notation, am is the Jacobi ampli- 
tude function, the solution of the sine-Gordon equation 
(3i| whose modulus is k. Plane waves are recovered for 
k — >• (H < H C 2), while the solitonic limit is reached at 
k ^ 1 4} (H > H cl ). For my = and mi,± ^ 0, the 
spin structure S^z exhibits zero magnetization, but the 
commensurate regions between soliton phase-slips cause 
the predicted NMR profile to exhibit a triple peak (see 
Fig. |4ji). This line profile is qualitatively the same as 
those which were claimed to indicate soliton formation 
in charge-density- wave materials 0, 14(1 . The correlation 
length £ m 4 is obtained via DMRG [37J, while k is close 
to one. To take into account the observed line shift of 




-0.2 0.2 0.4 0.6 -0.2 0.2 0.4 0.6 

Frequency [MHz] Frequency [MHz] n.02 

H > H 






-0.2 -0.1 0.1 0.2 0.3 0.4 0.5 

Frequency [MHz] 



0.01 



20 40 60 80 100 120 
„ Site along the leg 

||aBiCu 2 P0 6 -25T 
-Q — \ = 4, k = 0.9 
■S <S Z >=0.014 

- 

-0.2 -0.1 0.1 0.2 0.3 0.4 0.5 
Frequency [MHz] 



e) 



FIG. 4: (Color online) a) Simulated /(/). Left: with 77, e = 
and qi = 0.5 ±<5 (S ~ 0.04 and q is in reciprocal lattice units). 
Right: with e — 0, rj = 45° and q = qi (blue) or q — 2qi (red). 
m± was taken such that the widths are comparable with the 
data in panel e. b) Points are S z (i) values from DMRG, 
while solid lines are fits to the free particle state 5f ree ;(ip) 
{Slot — 8) an d soliton S^ ol i(ip) state (St ot = 1). c) Top view 
of the ladder in the C-phase predicted via DMRG calculations 
for H > H C 2 and D = 0. d) Expected NMR lineshape for a 
soliton lattice with mo = and mi,± 7^ and 77, e — 0. e) 
Soliton lattice as in d), but with mean moment (S z ) ^ and 
B Lp = 0.5. 



the NMR spectra, we have to consider also the longi- 
tudinal magnetization. We take mi t \\(ip) = S^ oll (ip), 
replacing Sf Teel (ip) with the periodic Jacobi function 
Ai t pdn[ip/(k^),k], which is a constant for k — > [3^ |. 
The resulting profile, as shown in Fig. mimics satis- 
factorily the experimental data. 

Based on the comparison between the observed and 
the calculated NMR lines, our interpretation of the IC 
phase applies to the field-induced order for small mag- 
netic moments at H c \ < H < H C 2 and H || b. When 
the frustration-driven short-range dimer order is broken, 
a lattice of QDWs appears. Our model also predicts a 
second-order transition at H C 2 3J| which, most likely, 
corresponds to the one observed in BiCt^POe at H2 — 
35 T using the same field orientation Based on the 
preceding discussion, the IC-C character of the transition 
can be verified via future NMR experiments. For H if( b, 
the g-tensor and the rung-DM anisotropy are more in- 
fluential and therefore will have to be included to model 
the complexity of the overall phase diagram [l8l ]. 

In summary, the field-induced order of the frustrated 
zigzag ladder Bi(Cui_j;Zn x )2P06 has been investigated 



via 31 P NMR with H \\b. NMR data and DMRG calcula- 
tions are combined to describe the magnetization process 
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as the formation of a lattice of QDWs. For x > the per- 
sistence of the ordered phase, an enhanced critical field 
and enhanced low-temperature spin fluctuations in the 
low-held region are reported. 
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SUPPLEMENTARY MATERIAL 

DMRG study of the soliton structure and of 
commensurate-to-incommensurate transition in 
BiCu 2 P0 6 

This is a short account of a DMRG numerical study of 
the ladder Hamiltonian introduced in Ref. It serves 
to complement the main results presented in our work. 
Here we focus on (i) the nature of the C-IC transition, 
and (ii) the structure of the correlations emerging as the 
spin system is partially polarized (at H > H c {), with the 
ultimate goal of presenting a phenomenological picture 
that can be compared to experiments. All calculations 
(unless specified) were performed on two-leg spin ladders 
of up to 128 rungs, and using 600 DMRG states. 

We start by recalling that when J 2 = 0, the geometry 
of the system shown in Fig. 1 of the main text, resem- 
bles that of two antiferromagnetically coupled sawtooth 
chains. Therefore, the physics of sawtooth chains, de- 
scribed in detail in Ref. [2|, may well be closely linked 
to that of BiCu2P06- For example, the classical ground 
state expected in BiCu 2 POe for J 2 = J 2 /2 has an in- 
commensurate spin order characterized by a wave vector 
q c = (0.608,1) close to 2/3, the value of coupled saw- 
tooth chains [2j. 

At a quantum level-the sawtooth chain exhibits spon- 
taneous dimerization [2|. However, when the spin chains 
are coupled in the transverse direction to form a ladder, 
the rung interactions can dominate 0], and the ground 
state of the system is better described as being formed 
by strong rung singlets. A strong rung-coupling regime 
has also been proposed for BiCu 2 P06 

To study the system correlations we define the dimer- 
ization operators along the leg and rung directions as 
Ao g (A,r) = S x , r ■ Sa, ? . + i, and D rung (r) = Si. r • S 2 ,r, 
respectively, with A = 1,2 being the leg index. The 
dimer order parameter is then defined as Odimcr(^ = i) = 
(D(i)) — (D(i — 1)), with i the site-index along the leg. 
From our calculations, both (D lung (r)) and (D\ cg (r)) are 
constant (except near the open boundaries), hence ruling 
out the presence of a spontaneous dimerization. In the 
unpolarizcd (H < H c {) case, we find that (Z? rU ng) is much 
stronger than {D\ Bg ) (not shown), in agreement with Ref. 
[H- One can go one step further and study the dimer- 
dimer correlator C(i,j) = (D(i)D(j)} - (D(i))(D(j)). 
Fig. UK shows rapidly decaying rung-rung correlations 
(calculated from the chain center). This fast decay is con- 
sistent with a state of weakly correlated rung singlets. As 
also shown in Fig. [5^, the dimcr-dimcr correlations along 
the leg direction indicate a cloud of sharply defined, but 
rapidly decaying staggered oscillations in a columnar pat- 
tern 0| . 

To break the rung singlets and create a nonzero magne- 
tization, a critical magnetic field H c \ is required. In this 
case, as shown in Fig. 4b of the main text, solitons are 
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FIG. 5: a) Dimer-dimer correlations in the ground state of the 
BiCu2POe spin-ladder model described in the text, calculated 
relative to the center of the chain. The rung-rung correlations 
and the dimer-dimer correlations, both along the same leg and 
between opposite legs are shown, b) Dimerization pattern 
along the legs and rungs of a polarized ladder with a pinning 
field at the center of one leg. c) Spin texture surrounding 
the soliton structure. The cartoon in the inset illustrates the 
classical spin texture and the structure of the spin correlations 
(thick lines indicate stronger bonds). 



formed and an IC S z modulation is present. Intuitively, 
solitons can be visualized as spin-i/2 QDWs breaking the 
short-range leg dimer order present for H < H c \ and vis- 
ible in Fig. [5^. Soliton pairs will then behave as triplet 
excitations with two "up" spins on the same rung. These 
pairs are weakly correlated along the leg direction and 
propagate coherently as hard-core bosons (see also Ref. 
[lj). Solitons in BiCu 2 P06 behave as weakly interacting 
extended objects, as shown in Figure 4b of the main text. 

However, as we shall see below, the results from our 
DMRG simulations indicate that the detailed spatial 
structure of the domain walls is more complex. A care- 
ful inspection shows the role of the incommensurate spin 
correlations with q rj 2/3, which are present as a quasi- 
classical spin texture superimposed on top of a dimer- 
ization pattern. To better characterize numerically the 
spatial profile of solitons, we added a pinning field at 
the center of the chain. The pinning field localizes a spin 
pointing in the "up" direction. Correspondingly, the spin 
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FIG. 6: Spin structure factor for a ladder with total spin S z — (a) and S~ = 2 (b). As magnetic field increases, a sharp peak 
develops in the transverse component, centered at q — q c . The inset in (c) magnifies the region around q c . Panel (d) shows 
the spatial structure of the correlations for S z = 10 as a function of the distance from the center of the chain. Notice that, in 
calculating the structure factor, the (S z (i)){S z (j)) contribution has been subtracted. 



on the same rung but opposite leg spontaneously points 
in the "down" direction. The surrounding background of 
spins arranges itself in a classical structure with q = q c . 
The excess spin "up" in the neighboring chain is not lo- 
calized on the same rung, but it spreads along the second 
chain over the width of the soliton, preferentially on the 
two sites to the left and right of the center, at an angle 
of <~ 120° obtained from the value of q c (see inset of Fig. 

Et). 

The soliton structure emerging from our calculations is 
unusual, and differs from those present in systems with a 
commensurate spin order. Moreover, due to the different 
couplings (J2 7^ J2), the profile of the soliton depends on 
whether it is situated on even or on odd sites of the leg, 
similarly to what happens in the sawtooth chains. No- 
tice that the reference position in our calculations corre- 
sponds to an odd site on one of the ladder legs. Results 
for even sites, although quantitatively different, are qual- 
itatively equivalent and hence are not shown here. 

In Fig. [6] we characterize the spin-spin correlations 
accompanying the FI phases. In the S z = case, the 
system is quantum disordered with a Lorentzian-like spin 
structure factor S(q) in reciprocal space (Fig. [6^). Once 
the critical field H c \ is reached, S z > and the ID system 
starts to develop quasi-long range order. The latter is 
manifest in the calculations with a sharp cusp in the spin 
structure factor, centered at q = q c (see Fig. [(J) and[Bfc). 
We also observe in Fig. IBJi that the correlation functions 




site i 



FIG. 7: Vector chiral correlations along a leg of the ladder, 
measured relative to the center, for a system of length L — 64. 
The inset shows the magnitude of the correlations 10 sites 
before the end of the chain, for different magnetizations. 



decay more rapidly for the spin components longitudinal 
to the applied field than for those perpendicular to it. 
The latter is a common situation in the gapless phase 
of two-leg ladders [1] and implies that in the presence 
of interladder interactions J', magnetic order will first 
involve the perpendicular component [f|. It has been 
shown 0] that long-range vector chiral order (breaking a 
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discrete parity symmetry Zo) can develop even in ID. A 
possible instability towards helical order in case of J' 7^ 
can be identified by numerically computing the chiral 
correlations (k(i)k(j)) with the chiral operator ft, = (S^ x 
S;+i) z . In Fig. [7] we see that vector chiral correlations 
are short-ranged in the region H c i < H < H C 2 , whereas 
a more classical cone is the ordered state realized in the 
3D case for H > H C 2. 

For H c i < H < H C 2, an incommensurate quasi- 
classical spin order develops around each QDW as the 
density of solitons increases. As in the case of CuGcC>3 
0], the interactions among solitons are supposed to af- 
fect mainly the soliton-soliton distance (and, in turn, my) 
but not the correlation length £. Therefore, at a critical 
density to|j corresponding to the inverse soliton width, 
namely m?, m l/2£ = 1/8, where £ « 4 is the correla- 
tions length, the solitons start to overlap and, due to the 
soliton-soliton repulsion, become more localized. In our 
system, we find that this corresponds to the development 
of long range chiral correlations at H C 2 . 

Due to simple energy considerations, the longitudi- 
nal magnetization parallel to the field develops predomi- 
nantly on one of the sublattices (even sites on leg 1 and 
odd sites on leg 2 — see Fig. 4b in the main text). The 
longitudinal magnetization is visible in the spin-spin cor- 
relations as a small commensurate feature with q = (1,1). 
It is not seen, though, in the structure factor, since the 



(S z (i))(S z (j)) contribution has been subtracted. Notice 
that this term would also have lead to a peak propor- 
tional to Sf ot at q = 0. 

To bear out the features associated with the soliton 
structure, we have ignored the effects of DM interactions 
in the above discussion. Adding them would add another 
degree of complexity to the problem Q. 
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